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INTRODUCTION 


This annual report outlines the major accomplishments of the Center 
for Remote Sensing and Cartography (CRSC) since the annual report was 
submitted In October 1982, with reference to the semiannual report of 
March 1983. The past year has been characterized by Important progress 
In the research of digital processing techniques, completion of projects 
Involving Integrated remote sensing and environmental analysis, and 
Interesting developments for new and continuing projects. Our computing 
strength is stronger than ever; we continue to rely primarily on the 
University of Utah Research Institute's Prime computer and NASA's Earth 
Resources Laboratory Applications Software (ELAS). Me continue to maintain 
and develop the software in-house, and through interaction with the newly 
formed ELAS Users Group. 

The staff at CRSC has gone through some changes over the last year. 
John Merola has succeeded Richard Jaynes as manager, and continues to 
serve the vital role as leader of the digital processing section. 

Additions to the staff Include two PhD candidates In Geography: Kevin 
Price, whose master's degree In Range Science adds substantially to our 
ecological studies, and Douglas Wheeler, who has a special Interest In 
geographic Information systems. Four work study students from the 
University of Utah have also been added; Steve Brower, from Mining Engin- 
eering; Jeanne Pearson, with a master's degree In Wildlife Science; Paul 
McDonald, a senior In Computer Science; and Gordon Douglass, with an 
Interest In computer mapping. Saundra Buckley continues to serve as 
secretary of the center. Wc* feel that we have a well-rounded staff with 
the ability to accomplish the goals we have set before us. 


CRSC FOCUS 


CRSC continues Its applied research In dealing with an array of 
environmental and resource problems In the I nte mountain Region. Our 
effort continues to focus on technical Innovations, as guided by the 
specific needs defined In each project. More and more, we are conscious 
of the value of dealing with each project as a team effort between CRSC 
specialists and the resource specialists. In this way, we build a bridge 
between CRSC and the resource manager, which adds credibility and useful- 
ness to our results. 

At the root of each project, the driving force Is the technical 
growth we can add through research developments. This approach leads to 
a continual CRSC concentration on alpine environments and arld/semlarld 
environments, the dominant landscapes of the region. This means that we 
must continue to be attuned to a variety of environmental systems. Includ- 
ing vegetation, soil, hydrology, geomorphology, and agriculture. In every 
project, the Interrelation of the various systems: Is not only evident, 
but significant In the problem-solving approach. This approach also leads 
to a continual stimulation to build strength In ancillary digital data, 
digital terrain data, and geographic Information systems (CIS) developments. 
Integrating them Into remote sensing techniques in a variety of ways. 

CRSC continues to pursue a course of merging remote sensing with CIS 
analysis. Ancillary data are vital In their own right as a decision-making 
tool, and also as a key to refinement In the classification of digital 
remote sensing data. Thus, we are moving toward an Integrated procedure. 
Invoking digital ancillary data and remote sensing data at appropriate 
levels, adding to the technical strength of CRSC. 




Folded Into this Integrated procedure Is the necessary field Investi- 
gation, photo Interpretation, map and archival data analysis to provide 
the resource manager with the comprehensive documentation necessary for 
sound decisions. 

We do not Intend to usurp, but to augment, the professional manager's 
role, filling those steps he Is unable to fill otherwise. Our staff 1$ 
well trained In field, laboratory, and library procedures In resource 
analysis as well as In remote sensing per se. Recently, we have been 
called on to perform such comprehensive evaluations, e.g., the Parker 
Mountain study for the Utah Division of State Lands and Forestry, and the 
Hueboldt River study for Nevada Fish and Game. In these broad studies, as 
well as In projects more restricted to techniques of remote sensing, CRSC 
Is ready to bring In the agency professional to the fullest extent. We 
wish to have agency personnel fully aware of the capabilities and 
limitations of the role of remote sensing and GIS analysis, so that their 
professional judgement Is not abdicated to our laboratory or to the 
"machine." Meanwhile, CRSC adds another Increment of basic as well as 
applied technical capability. 
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RECENT FUNDAMENTAL RESEARCH DEVELOPMENTS 


The research developments presented here are those primarily developed, 
or concluded, over the last year. We begin with two research papers, produced 
out of project applications, and proceed to describe developments In use. 

Identifying Ecological Range Sites Through Remote Sensing 

This procedure was developed during the applied project on Parker 
Mountain. The project Involved Inventory and analysis of the rangeland; 

This procedure Involves the use of map Information derived from both high 
altitude color Infrared photography and Landsat digital data. The maps are 
an Integration of the two remote sensing sources along with soils, geology, 
and precipitation. Identification of the ecological range site for a given 
area requires an evaluation of all relevant environmental factors. These 
environmental factors combine to give a site the potential to produce char- 
acteristic types and amounts of vegetation. The utilization of the photo- 
graphy and Landsat digital data, coupled with digital ancillary data on soils 
geology and precl citation. Is shown to be a means of achieving greater 
efficiency and accuracy In Identifying ecological range sites. This paper Is 
attached as Appendix A. 

Contextual Analysis of Landsat Spectral Signatures 

This paper (an abstract of which Is included as Appendix B) Is the 
conclusion of several years of development and testing. The procedure 
presented has proven Itself In many projects and has become a "standard" 
approach to the classification of digital Image data at CRSC. 

The procedure consists of a sequence of three multivariate statistical 


routines to aid In detecting similarities and differences In spectral 


signatures. A wore precise computer classified land cover wap Is produced 
by using and understanding additional spectral classes that are represented by 
the signature. The sequence of statistical routines used are: principle 
components analysis, cluster analysis, and discriminant analysis. 

Principle components analysis Is used to transform the data set 
(signatures) Into Independent variables to remove colllnearity. The 
variables now meet the assumptions of the statistical routines that follow. 

The output data are factor scores for two components representing the 
original data set. The factor scores are then used In a cluster analysis to 
group spectral signatures according to a similarity Index. Finally, the 
factor scores and group clusters are used In a discriminant analysis. A 
two-dimensional scatter plot produced In the discriminant analysis allows 
one to obtain a graphical view of the spectral context In which a particular 
signature Is found. 

The use of discriminant analysis. In combination with examination of 
spectral signature plots and ground Information, has been a key element In 
achieving good results from the unsupervised approach to Landsat data 
analysis. 

Spectral and Spatial Contextual Analysis of Digital Imagery 

At CRSC, we have developed and put Into operation a technique for acquir- 
ing Information about the context of a spectral signature. The procedure Is 
analogous to analyzing a word In a sentence; take a word out of a sentence and 
It can be hard to understand the meaning of that word unless It Is related 
to the other words In the sentence. Our procedure relates each spectral 
signature to all others acquired In a particular study area. The technique 


provides Information about the Interrelationships, specifically the spectral 
context of each signature In the set of signatures. Now that we are able 
to gather spectral contextual Information about the signatures used to 
classify digital Image data, we decided the next logical step was to Include 
spatial contextual Information for classification Improvement. We believe 
this technique will be most useful for high resolution digital Image data, 
such as Thematic Mapper (TM) data from the Landsat 0 series. Currently, the 
computer program has basically been developed and Is operational. Testing 
has begun on a Thematic Mapper Simulated (TMS) data for Salt Lake County, 
1982, but conclusive results are as yet unavailable. 

Processing of Thematic Mapper Data 

Several techniques have been evaluated to date on TMS data covering 
Salt Lake County. The purpose of this evaluation Is first to prepare CR5C 
for the use of Landsat 4 TM data, and secondly, to compare TMS data to MSS 
In the urban environment of Salt Lake County. Our evaluation has focused 
on the spectral bands most useful for classification of land cover In the 
stu4y area. The results show that bands 2, 3, and 4, In combination with 
band 7, are the most accurate for the broad application of general land 
cover mapping. 

Evaluation of Change Detection Routines 

We have committed much effort to the detection of agrlcuiture-to-urban 
land use change using multi temporal MSS digital data. Although Landsat MSS 
data Is a relatively coarse tool for discriminating categories of change In 
the urban size plots. Its availability over the 11 years of the Landsat 
program prompts a thorough test of Its power to detect change. A recently 
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presented research paper Is attached as Appendix C. Many techniques have 
been tested and results are nearly completed. The results will be 
published In the upcoming paper: "Evaluating Landsat MSS Digital Data for 
Change Detection; on the Urban Fringe." 

Developments In Computer Programming 

In the process of accomplishing our goals In the applied research 
projects we are Involved In, new computer programs are being developed. They 
are being developed to augment and become part of our software package, ELAS. 
The following Is a short list of programs developed at CRSC: 

1. PLYI - polygon Insert. This module allows the Insert of polygons 
Into an ELAS data file with no effect on the data outside the 
polygon. 

2. PCHG - pixel change detection. This module Is used to find changes 
In pixels from one date to another. It Is used for continuous 
data and Is specifically designed for use with albedo change 
detection. 

3. PMIN - pixel minimum. Find the minimum pixel value for any or all 
channels In an ELAS data file. 

4. CNTX - contextual reclasslfler. This program reclassifies pixels 
based on spectral and spatial context Information described 
previously. 

Fully documented copies of these programs are, or will be, available to any 
ELAS user Interested In them. 

Another program, which Is now part of ELAS, that has received some 
development but deserves special consideration Is ACTS (Accuracy of Classi- 
fication Table Builder)* The reason ACTB deserves special consideration Is 
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that It hat become the pivot cf a “standardized procedure" for the calibration 
and verification of naps which result from the Investigation Into the pro- 
cessing of digital Image data here at CRSC, At stated In the ELAS manual 
(Report No. 183) User's Guide* Section S: 

"5.16.5 ACTS - Accuracy of Classification Table 

ACTB compares the results of a classification with 'ground truth' 
or verification data. ACTB outputs a table that shows class frequencies, 
percentages, percent correct, omission errors, t*nd commission errors as 
a result of the comparison between the verification data and the 
classified data." 

At CRSC, our "standardized procedure” for calibration begins by select- 
ing, from air photos, ground verification plots. These plots are then 
digitized and are compared to the several digital Image processing techniques 
being investigated. This Is an Iterating process in which the most accurate 
technique 1$ selected. 

ACTB Is also used to verify the final map accuracy by adding to, or 
using separately, a set of new verification plots. These new plots are used 
to make a statement of the overall accuracy of the final maps. This has 
proved to be an Important step In bridging the gap between CRSC and the 
resource manager. By using an accepted standard approach In stating map 
accuracy, we have found coranon ground In which to work together with the 
resource managers. This is as important to our research strategy as It Is 
to users of our work. 
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APPLIED RESEARCH PROJECTS COMPLETED 


Three projects have been completed since our *mua1 report was prepared 
In October 1982. One of these projects has been coiRpleted since the semi* 
annual report In March 1983. Highlights of each project are presented below. 

Salt Lake County Land Use 

Because of rapid urban expansion along the Wasatch Front, many Impacts 
are created, l.e., loss of farmland, changing water needs, transportation 
corridor demands, etc. Accurate update maps of urban expansion and land use 
change are needed by many agencies. This project Is an attempt to answer to 
that need. It consists 0 / two parts: (1) land use Inventory, current to 
1982, using color Infrared (CIR) photography, establishing a sound baseline 
for change detection, and 1 2} an experiment to evaluate the effectiveness of 
Landsat MSS digital data to detect, measure, and map change In this complex 
environment. The CIR Inventory Is complete. Five copies of the report were 
forwarded to NASA Headquarters on October 5, 1983. Twenty- two quadrangle 
overlays at 1:24,000 have been presented to the Utah Division of Water 
Resources. The maps have been carefully checked, drafted, and turned over 
to the Division. 

The new overlays are currently being used In their own right by the 
Division for Water Allocation Modeling. Other departments of state and 
local government are accessing them for their own purposes. For CRSC, the 
completed maps serve as the base for phase two of the study: evaluation 
of digital techniques for rural -urban charge detection. Several technical 
Innovations are being tested by CRSC In a project underway to detect change 
on the urban fringe. 
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Parker Mountain State Land Block 

To test the effectiveness of Landsat data In detecting pasture quality 
and potential for grazing production. CRSC completed a comprehensive 
analysis of 39.292 acres of state land In south-central Utah. Five copies 
of the report were forwarded to NASA In March 1983. 

This project established another milestone for CRSC. by Integrating 
digitized ancillary geological data and elevation data with Landsat digital 
data, using ELAS program routines. Accuracies of classification of range- 
land Increased significantly. Using Landsat MSS data alone, accuracy of 
classification was 56X. By digitizing various basalt units taken from the 
state geologic map. and critical elevation lines as surrogates for precipi- 
tation variation, accuracy of classification Increased to 77X, a substantial 
21% over the best solution from spectral data alone. 

The technical advancement to CRSC from this project was significant. 

In addition to adding ancillary digital capabilities, we developed Improved 
sampling strategies for accuracy, verification, and gained further experi- 
ence in principle components and other techniques. 

Aspen-Conifer Succession Study 

In a joint project with the Intermountain Forest and Range Experiment 
Station (U.S. Forest Service), a quantitative index was developed which 
correlates with apparently stable and serai aspen forests. Several classes 
of stable/seral aspen forests have been mapped in the Bear River Range of 
Utah and Idaho. Since aspen canopies tend to obscure understory conifers 
for early serai forests, a second date, using Landsat data taken when aspen 
trees are leafless, has been used. This study explored the extent to which 
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a two- date Landsat MSS analysis Improves the detection of aspen/conifer 
forest nixes. By using two dates, our overall accuracy Improved only by 
4%, from 68% to 72%, but the accuracy of stable aspen Improved from 79% to 
90%; and late serai Improved from 41% to 65%; conifer remained at 85%. As 
these are very Important classes, from a forest manager's point of view, we 
feel our goals for this study were met, and then some. 
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APPLIED RESEARCH PROJECTS UNDERMAY 


At this time, there are eight projects underway. These projects vary 
In their dimensions from small* specific Investigation, fine tuning our 
recent developments* to the larger multifaceted Investigation* encompassing 
many techniques and requiring the development of new capabilities. The 
Importance of these applied research projects cannot be over-emphasized. 

CRSC could never have accomplished as much as we have In remote sensing 
research without the ability to work closely with, and provide results to, 
the user community, l.e., resource managers. Highlights of each project are 
presented below. 

Evaluating Landsat MSS Digital Data for Change Detection on the Urban Fringe 
This Is a continuation of the Salt Lake County study In rural-urban 
change detection. The objective Is to compare several digital procedures 
such as: band five differencing, data transformation using the Kauth trans- 
formation, and other algorithms for detecting land conversion patterns. 
Categories of land use change are selected on the basis of water consumption. 
The question is, how well can MSS data Identify these categories and patterns 
of change? Several algorithms are being compared, utilizing both raw data 
and preprocessed data. Verification of results Involve high quality color 
Infrared (CIR) photography and field observation. The dates involved in 
the test are 1972, 1979, and 1981. "Change Maps" are produced at 1:24,000 
scale. This digital analysis should help develop a procedure for wider 
area application and for subsequent updating without the necessity of 
repeated CIR photography. 
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Thematic Mapper Simulator Data Analysis 

The thematic mapper simulator (TN5) Is a multi spectral scanner flown 
from an aircraft to simulate the TN data that will be available from 
Landsat 4. The goal In this analysis 1$ to compare MSS data from earlier 
Landsat satellites and IMS data to Identify the possible role of TM and Its 
use In the future. 

It has been stated by several authors that an Improved spatial and 
spectral sensor, such as the TMS sensor, does not necessarily Improve 
classification accuracies when processed the same as MSS data. The pro- 
blem, as stated by Warton (1982), Is as follows: 

"To Illustrate the problem, consider the task of distinguishing 
between commercial and residential land-use classes in 
[spatially] high-resolution remotely sensed data. The straight- 
forward per-pixel spectral classification of the data can Identify 
only the spectrally dissimilar ground cover classes such as pave- 
ment, lawn, trees, roof. These components are common to both the 
commercial and residential land-use classes, making It difficult 
to distinguish between them on a per-pixel basis. It has been 
suggested by several authors that the solution to this problem 
is the use of contextual information (i.e., the local frequency 
distribution of components surrounding each pixel)." 

A contextual reclassifier has been developed to address this problem. 
Spatial and spectral contextual Information is applied in the procedure. 

The spectral contextual information is supplied by the statistical analysis 
of the signatures, as described previously. The two-dimensional scatter 
plot produced in the discriminant analysis is based on discriminant scores 
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derived from the analysis. These discriminant scores supply the spectral 
contextual Information. The spatial Information Is supplied by a moving 
window of specified dimensions* and a similarity matrix (derived from the 
discriminant scores) which evaluates the center pixel of the moving window. 

The center pixel Is then reclassified If the specified parameters are met. 

This procedure will be tested on MSS as well as TM$ data. 

Humboldt River Riparian Habitat 

The Nevada Department of Wildlife recently contacted CRSC regarding the 
need to obtain accurate, up-to-date maps of riparian habitat along a 
segment of the Humboldt River In north-central Nevada. Wildlife habitat Is 
significantly affected by various activities along the river, 1.e. t channelization 
and spraying or burning of willow areas. The primary objective of the study 
1$ to prepare 1:24,000 scale topographic map overlays of riparian habitat 
and other significant features for a stretch of river which flows through 
portions of land represented on eight U.S.G.S. quadrangles. This mapping 
project will explore the extent to which environmental changes due to 
management activities may be monitored over time. 

The primary purpose of the study is to provide a spatially-oriented 
Information base for wildlife population and habitat management for this 
particular riparian zone. This Information base will be a key element In 
the development of Department policies and public awareness regarding 
proposals and actions for river channelization and phreatophyte 
control. The study will also demonstrate the limits of using available 
aerial photography for obtaining necessary riparian habitat Information. 

Due to unusual amounts of precipitation last winter, the Humboldt 
River riparian area was flooded to a greater depth and for a longer than 


- 14 - 



usual period this spring. This caused somewhat of a delay In getting onto 
the ground and beginning preliminary field studies. This rare event pro- 
vided the opportunity to record the water Inundated areas along the river. 
On June 8* most of the Humboldt River was photographed, from Deeth to 
approximately IS miles east of Battle Mountain. The photographs were taken 
from 3.000 feet above ground level using a panoramic camera (Envlropod) and 
color Infrared film. Please refer to the semiannual report, March 1983, 
Exhibit D, for further detail about this study. 


Agriculture and Wetlands Inventory In the Lower Sevier Basin 

The objective of this study Is to classify and map agriculture-related 
land cover and wetlands of the Lower Sevier River Basin (l.e., below Piute 
Reservoir). The mapping will be accomplished primarily by Interpreting 
available high altitude color Infrared (CIR) photography. Map overlays 
will be prepared to match existing U.S.6.S. topographic quadrangles. The 
acreage for agricultural lands and wetlands will also be tabulated. 

This study will provide basic Inventory information about water-related 
land use for application In basin-wide water budgeting by the Division of 
Water Resources. The coordinated mapping of wetlands and agricultural lands 
will also assist farmers and federal and state agency personnel In better 
understanding the relationship between wetlands and farms. Such Information 
will promote more Informed farm and ranch planning and assist agencies such 
as the S.C.S. in carrying out their responsibilities under legislation and 
regulations which protect wetlands. Finally, this study will provide a 
basis for documenting changes In agricultural land use since the last 
inventory, some 20 years ago, and will permit greater detail In future 
monitoring of resource uses In the Sevier River Basin. 



EPA Envlropod 

CRSC has signed an agreement with EPA and state agencies to test and 
evaluate the role of small format panoramic camera use In environmental 
monitoring and assessment. CRSC will serve as the technical advisor. 
Involvement In this program Is Ideal for CRSC since each Envlropod appli- 
cation could provide support Information for our Landsat data and allow 
for evaluating the comparative utility of various stages of remote sensing 
data. 

The Envlropod Is a compact, self-contained, panoramic, two-camera 
reconnaissance system designed to be secured to widely available aircraft 
(Cessna 172 and 182). We have flown a total of 20 missions to date. The 
majority of applications have related to flooding, but land cover/land use 
Inventories have also been flown. These land cover/land use Inventories 
have given us the necessary experience to use the Envlropod to support 
projects Investigating the use of digital landsat data. The 
comparative rangeland study In Rush Valley, Utah will be our first 
investigation and application of this approach. Another Investigation of 
this approach and of other techniques will be performed In a study of 
desertification In west-central Utah. A summary of these applied research 
projects follows. 

Comparative Condition Study of Desert Rangelands In Rush Valley, Utah 

One of the key Issues in rangeland management is the uncertainty of 
stocking capacity on specific pastures. Closely related Is the changing 
nature of that capacity through time. Range managers have attempted to 
resolve the problem by “measuring" condition and trend of the range. One 
of the difficulties Is that the transect site (and the specific plants) 
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obtained through NASA's Landsat satalllta lyitn from tha 
EROS Data Cantar. For data processing, tha Center is 
utilising tha "ELAS" package of computer software routines, 
developed by NASA's Earth Resources Laboratory, which is 
operational on tha Univarsity of Utah Research Instituta's 
PRIME coaputar . 

Tha priaary ratlonala for parfonsing digital process- 
ing of MSS data begins with tha assumption that different 
typas of ground covar hava diffarant patterns of raf lection. 
It is also assumed that these spectral patterns are 
sufficiently unique to naka diffarant ground cover typas 
consistently distinguishable froa one another using 
statistical classification techniques (Hutchinson 1982) . 


Standard Digital Data Analysis 


Ground cover characteristics nay be analysed using 
digital laage data from a variety of sources: Landsat 

multispectral scanner (MSS) ; Landsat thanatle nappcr (IK) ; 
airborne thematic mapper simulator (IMS) ; NOAA satalllta 
sensor, advanced vary high resolution radiometer (AVHRR); 
etc. The description of methods in this section applies to 
the analysis of digital image data from any such source. 


Any digital image data must first be reformatted to 
make them compatible with processing hardware. Next, a 
grey level map is produced which encompasses the study area 
and some of the surrounding area. The size of the grey 
level map is determined by the number of features present 
to be used as ground control. As water bodies typically 
provide the sharpest and most spectrally distinctive 
features for control, a near infrared band is generally 
used to produce the map. (See Stage 1 of Figure 1.) These 
ground control points are used later to geographically 
reference the digital data to remove the effects of earth 
curvature, spin, etc. (See Stage 3 of Figure 1.) At this 
point of data processing, the raw data may be modified by 
other analytical steps such as filtering, principal compon- 
ents analysis, etc. Thereafter, a program called "SEARCH" 
is utilized to generate statistics which characterize 
pixel groups having similar spectral features across the 
bands.. (See Stage 2 of Figure 1.) SEARCH is -a routine 
which is used to provide training statistics for a program 
called "MAXL," which classifies individual pixels into a 
class based upon each pixel's highest statistical probabil- 
ity of belonging to a given decs. The training statistics 
are derived from blocks or windows of data which correspond 
to areas on the ground containing the cover types of 
Interest. 


Once the study windows are selected, the SEARCH program 
examines each continguous six scan line (Landsat pixel 
matrix "row") by six element (pixel matrix "column") block; 
if the spectral data within the six by six block are too 
heterogeneous , the program will switch to the use of a three 
by three block of pixels. The statistics generated by 
SEARCH include mean pixel light reflectance values for each 
of the four bands, a covariance matrix, and a priori values. 


that salinization can be detected on Landsat based on an Increase In 
brightness of spectral response due to the salt content In the upper crusts 
of the soils. Similarly, erosion should also exhibit a brighter spectral 
response due to lack of vegetation in eroded areas. This knowledge will 
aid understanding the dynamics of desertification and thereby improve 
resource management of the environment. A second objective is to improve 
the interpretation techniques of satellite digital image data as it relates 
to arid environment studies by combining raw data spectral signatures with 
ancillary field data to derive a workable approach for delineating areas of 
salization and erosion. The resource-oriented capabilities of remote 
sensing, applications in computer mapping, the use of information in a 
multi-disciplinary Geographic Information System, and interrelationships of 
the physical and biological landscape will combine to provide significant 
information to address the problems of desertification. 

Mapping of Alpine Riparian Areas from Landsat MSS and Terrain Digital Data, 
Uinta Mountains, Utah 

A significant amount of investigation into mapping alpine riparian 
areas using Landsat digital data exclusively has already been accomplished. 

In this investigation, it was discovered that stands of coniferous trees 
that are associated with water, i.e., standing in very wet soil and classified 
riparian, were separated consistently from dry conifer stands. During the 
verification of the maps, it was also discovered that shadowed slopes were 
consistently being confused with riparian conifer stands. To resolve this 
problem, it was proposed that digital terrain data be used in conjunction 
with Landsat digital data to resolve the confusion. Due to the interest 
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of the U.S. Forest Service In the use of digital terrain data, they agreed 
to reduce the size of the mapping project In order to include the analysis 
of digital terrain data. 

At this time, we are still trying to Implement the software that will 
allow us to use c' gltal terrain data. We reported In our follow-on proposal 
and semiannual report that a modification to the Prime computer operating 
system had given us this capability. In testing the modification, it was 
discovered that It doesn't work. An alternative Is to acquire from the 
Forest Service the software necessary to process the data on the University 
of Utah's Univac 1100/61 computer. This software has been received and 
will be Implemented soon. The software called MULDEM will allow us to 
preprocess digital terrain data, on the less expensive Univac, and then 
bring the data (In a ready to use form) to the Prime. We will have to write 
a short program to allow ELAS to use this data. It is anticipated that we 
can finish this project by spring of 1984. 
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PROPOSED FUTURE RESEARCH 


Fundamental Research 

Software development Is an area where a major effort will be made to 
improve our current capabilities. The Improvement in software will become 
part of ELAS and will add to Its already extensive capability. 

On the topic of digital geographic information, which is vital to the 
future use of digital remote sensing data, we propose several areas of 
Investigation such as the ability to produce polygons from the pixels of a 
classified digital image map. This is important in developing the capability to 
access the recent developments in polygon manipulation and analysis of 
digital geographic Information. Specifically, we want to work closely with 
the state of Utah's Automated Geographic Reference (AGR) system which uses 
the Arc/Info software system of Environmental Sciences Research Institute's 
(ESRI) for polygon manipulation and analysis. At a recent meeting on 
October 21, 1983, of the Governor's Panel on Scientific Inputs to Land Use 
Planning, AGR Director, Riki Darling, disclosed the intent of the AGR to 
purchase a Prime computer and the ESRI Arc/Info software. Given our experi- 
ence In the research and analysis of digital geographic and remote sensing 
information to produce land cover/land use maps, the capabilities of ELAS, 
and our use of a Prime computer put us in a very good position to provide 
technical support to AGR. Discussions to this end have already begun 
between CRSC and AGR. 

Another area we intend to continue to investigate is the use of 
contextual information for the classification of land cover/land use from 
digital image data. We propose to develop the capability to create a 


previously undefined class of land cover/land use based on the contextual 
Information In a given window, or subscene, of a digital Image. 

Although CRSC and other facilities have been Investigating the use of 
digital map data in GIS analyses, there Is much still to be discovered. 

Once a map Is In the digital form, the digital data are then used In various 
types of GIS analyses. 

The ELAS program package allows us considerable flexibility In using 
the digitized maps for post-classification Improvement of MSS digital data. 
For example, our analysis of rangeland in the Parker Mountain study was 
greatly Improved using this technique. Polygons derived from Sv.il and 
geology maps were digitized for stratification of the classification. 
Selected contour lines digitized from topographic maps also added signifi- 
cantly to the classification accuracy. Combined, the digitized entry of 
soil, geology, and elevation data produced a classification accuracy which 
was over 21% more accurate than the maps prepared from spectral data alone. 

A new computer program was written to insert the polygons back into the map 
after changes were made. 

ELAS has also allowci us to use a very complicated GIS analysis for 
the combination of two dates of MSS data in the Aspen Succession study. The 
analysis was complicated in the sense that there were many possible changes 
that could have been made to any one pixel provided that it passed one or a 
combination of seven different tests. We are very pleased with the results 
obtained thus far. There are several other areas of fundamental research 
we intend to explore. For more information, please refe to our follow-on 
proposal submitted in April 1983. 
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PROPOSED PROJECTS FOR RESEARCH AND DEVELOPMENT 


Land Use/Land Cover Inventories and Data Integration to GIS 

In this project. It Is proposed that we develop the technique of 
assessing land use/land cover In a digital format in order to Integrate 
the data in the State of Utah's Automated Geographic Reference (AGR) 
system, l.e., GIS. A memo requesting proposals to accomplish this appears 
as Appendix D. Each of the three proposals is in a different geographic 
area within the northern part of the state. These areas have different 
land use/land cover types as well as various environments. We have chosen 
the North Cache study area for a first analysis and development of the 
techniques that will be used. 

Soil Loss Evaluation in Pinyon -Juniper Using a GIS Format in Utah 
We are exploring the possibilities of a project with the Soil 
Conservation Service to assess the capability of Landsat. digital data, 
coupled with various GIS inputs, to detect variations in pinyon* juniper 
cover that con al ate with erosion conditions. The experiment will be 
built around the Universal Soil Loss Equation. 


PROFESSIONAL MEETINGS ATTENDED BY CRSC PERSONNEL 


American Association of Geographers Rocky Mountain Regional Meetings, 
Boulder, Colorado, October 1983. The following papers were presented. 

Rldd, Merrill K. , "Hazards to Development on the Wasatch Front: An 
Application of Geographic Information Systems Analysis " 

Merola, John A., "Aspen Succession" 

ACSM-ASP Fall Convention, Salt Lake City, Utah, September 1983. The 
following papers were presented. 

Jaynes, Richard A., "An Integrated Remote Sensing Approach for 
Identifying Ecological Range Sites" 

Jaynes, Richard A., "Remote Sensing Evidence and Experts In the 
Courtroom" 

Merola, John A., Richard A. Jaynes, and Roy 0. Harnlss, "Detection of 
Aspen/Conifer Forest Mixes from Multi temporal Landsat Digital 
Data" 

Rldd, Merrill K., John A. Me~ola, and Richard A. Jaynes, "Detecting 
Agricultural to Urban Land Use Change from Multi temporal MSS 
Digital Data" 

Rldd, Merrill K. , "Hazards to Development on the Wasatch Front: An 

Application of Geographic Information systems Analysis" 

ELAS Users Group Meeting, Bay St. Louis, Mississippi, NASA/NSTL/ERL, June 
1983. John A. Merola made a presentation of research using ELAS. He was 
also appointed to chair the Documentation Committee. 

International Symposium on Remote Sensing of Environment, Ann Art.r, 
Michigan, May lv83. The following paper was presented. 

Merola, John A., Richard A. Jaynes, and Roy 0. Harniss, "Detection of 
Aspen/Conifer Forest Mixes from Multitemporal Landsat Digital 
Data" 

RNRF Symposium on the Application of Remote Sensing to Resource Management, 
Seattle, Washington, May 1983, Richard A. Jaynes attended. 

American Association of Geographers Annual Meeting, Den\er, Colorado, April 
1983. The following paper was presented. 

Ri»1d, Merrill K. , "Color Stretching and Digital Analysis of Landsat 
Data: An Inventory of Guayule in Mexico" 
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ASCto-ASP Annual Meeting, Washington, D.C., March 1983. The following paper 
was presented. 

Ridd, Merrill K. , "Preliminary Digital Classification of Grazing 
Resources in the Southern Chihuahuan Arid Zone of Mexico" 

Working session of Caribou National Forest, Idaho Fish and Game, etc., March 
1983. John Merola presented the aspen succession study. 

Meeting of Environmental Health Division, January 1983. Merrill Ridd was 
invited to present a summary of CRSC remote sensing activity. 

International Symposium on Remote Sensing of Environment Second Thematic 
Conference: Remote Sensing for Exploration Geology, Fort Worth, Texas, 
December 1982. Merrill Ridd attended. 

Utah meeting of ACSM-ASP. Richard Jaynes and Merrill Ridd presented a 
summary of CRSC projects. Held December 1982. 

Bear River Resource ?*nd Conservation District meeting in Logan, Utah, 
December 1982. Merrill Ridd was appointed to a committee to create a pilot 
project for comprehensive soil and water management planning, based on 
remote sensing and GIS. 

Utah Association of Geographers, November 1982. Merrill Ridd presented a 
paper on remote sensing studies of the Great Salt Lake. 

Merrill Ridd was appointed to ‘he Governor's Science Advisory Committee. 
Panel on Scientific Inputs to Land Use Planning 


Following is an updated list of projects and reports by CRSC under 
NASA (or partial NASA) funding. 



Exhibit A 


Table 1. CRSC Projects Supported in Whole or Part by NASA Granti NSG-7226 and NAGW-9S 


Project Short Title 

CRSC 

Report 

Agcncy/ies 

Agency 

Support 

Completed 

Project Impact 

PROJECTS COMPLETED 






North Ogden Hazards to 

78-1 

Weber County, 

Limited 

May 1978 

Adopted for Sensitive Ana Overlay Zone 

Urban Development 


No. Ogden City. 
Pleasant View 

in Kind 


Ordinance. 

Irrigation Detection by 

79-1 

Utah Div. of 

Limited 

April 1978 

Proven effective. Led to ground water study. 

Satellite, Iron Co. 


Water Rights 

in Kind 



Price River Bat.ii 

79-2 

U.S. Bureau of 

Limited 

May 1979 

Proven effective but demands real-time data 

Rangeland Response to 
Summer Rain 


Land Management 

in Kind 


for application. 

Korean Land Use II 

79-3 

Republic of Korea 

$102,000 

Aug. 1979 

NASA funding helped develop software; 
technical development. 

Snowpack/Runoff 

Exper- 

Utah Div. of 

Limited 

.. 

High correlations shown. Tabled until near 

Correlation 

imental 

Water Resources 
Soil Cons. Sent. 



real-time data is available. 

Guayula Inventory 1: 

80-1 

Mexican 

$90,000 

Jan. 1980 

NASA funding helped perfect computer 

Contrast Enhancement 


Government 



enhancement 

Uinta Basin Wetland/Land 

80-2 

Utah Div. of 

$10,000 

Dec. 1980 

Wetland management, weter allocation and 

Use 


Water Resources 
Soil Cons. Serv. 

$25,000 


management, agriculture resource planning. 

Snow Cover/Mule Deer 

80-3 

Utah Div. of 

Limited 

July 1980 

Landsat utility, tabled pending agency studies. 



Wildlife Res. 

in Kind 



MX Draft EIS Review 

81-1 

Utah Governor's 
Office 

ca. $2,000 

April 1981 

Used by Governor for MX policy and comment 

Guayula Inventory II: 

81-2 

Mexican 

$100,000 

Feb. 1981 

NASA funding helped develop classification 

Statistical Classification Routines 


Government 



routines. 

Farmington Bay Sh oreline 

81-3 

Utah Div. of 

In Kind 

April 1981 

Deterred proposed project which would have 



Wildlife Res. 
Great Salt 
Board 

Minimal 


damaged waterfowl habitat 

Irrigation Detection, 

81-4 

Utah Div. of 

Minimal 

May 1981 

Proven effective; led to Bear River Study; 

Iron County II 


Watar Rights 



leading to prosecution. 

Davis County Foothill 

81-5 

Davis County 

In Kind 

May 1981 

Being reviewed for edoption as the guideline 

Development 


Planning Comm. 
Four Corners 

$63,000 


for urban development control 



Regional Comm. 
EPA 208 

$18,000 





Weber-Oavis 
Several State 

In Kind 





Agencies 

Seven municipalities 

In Kind 
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Appendix A 


AN INTEGRATED REMOTE SENSING APPROACH 
FCR IDENTIFYING ECOLOGICAL RANGE SITES* 

Richard A. Jaynes 

Center for Remote Sensing and Cartography 
University of Utah Research Institute 
420 Chlpeta Way, Suite 190 
Salt Lake City, Utah 84108 

ABSTRACT 

A model approach for Identifying ecological range sites 
has been applied to high elevation sagebrush-dominated 
rangelands on Parker Mountain, in south-central Utah. The 
approach utilizes map Information derived from both high 
altitude color infrared photography and Landsat digital 
data, integrated with soils, geological, and precipitation 
maps. Identification of the ecological range site for a 
given area requires an evaluation of all relevant environ- 
mental factors which combine to give that site the poten- 
tial to produce characteristic types and amounts of 
vegetation. A table is presented which allows the user to 
determine ecological range site based upon an integrated 
use of the maps which were prepared in this study. The 
advantages of identifying ecological raDge sites through an 
integrated photo lnterpretatlon/Landsat analysis are 
discussed. 


INTRODUCTION 

The ecological range site concept has become the most 
widely used foundation for range management In recent years. 
An ecological range site (hereinafter "range site") is a 
"distinctive kind of rangeland that differs from other kinds 
of rangeland in its ability to produce a characteristic 
natural plant community" (S.C.S. 1976). The natural plant 
community or climax community is that assembledge of plants 
that v*ould eventually occupy a site in the absence of 
abnormal disturbances and physical site deterioration. 

Range sites are derived from the analysis of vegetation 
composition (by dry weight) from relict sites with similar 
soil, climatic, topographic, and geologic characteristics. 
Plant association tables are prepared and analyzed for 
significant differences in the kind of dominant species and 
species groups, proportionate make up of dominant species 
and species groups, and total annual production. Thus, 
definition of range site and designation of an area as being 
that range site, provides a single expression of all 
environmental factors responsible for the development of 
that range site. The range site concept has, therefore, 
become a key component in the use, development, and rehab- 
ilitation of rangelands. Although a given range site model 
may oversimplify the inherent variability in nature, it is 

* This study was supported by the Utah Division of State 
Lands and Forestry, and the National Aeronautics and Space 
Administration (Grant NAGW-95). 


nevertheless s velueble end edsptable tool for developing 
rengelend management plans. 

If rangelands were In relict condition, one would 
merely need to map vegetation types to identify range sites. 
Since many rangelands have received intensive use, range 
site identification is achieved through an integrated 
analysis of vegetation (considered a temporal attribute) and 
spatial attributes. 

The traditional approach to mapping range sites is 
field labor intensive, with reliance on large scale black 
and white aerial photography for vegetation and geomorphic 
analysis. Spatial attributes are also obtained from avail- 
able soils, geologic, topographic, and precipitation maps. 

The objective of this study was to explore the utilization 
of high altitude color infrared ("CIR*’) photography and 
Landsat digital data as a means of achieving greater 
efficiency and accuracy in the identification of ecological 
range sites. This study is part of a comprehensive study 
to analyze and map rangeland resources of the Parker Mountain 
study area (described below) for the Utah Division of State 
Lands and Forestry (Jaynes 1982) . 

STUDY AREA 

Efforts to identify range sites in this study were 
carried out on the Parker Mountain State Land Block in 
south-central Utah. The study area occupies over 45,000 
acres of high elevation (8,600-9,800 ft.) rangeland on the 
western edge of the Awapa Plateau, an eastward sloping 
plateau covered with various types of volcanic flows and 
deposits. The Parker Mountain study area is characterized 
by rolling hills covered with mountain big sagebrush and 
islands of aspen forests on its western half, and, on the 
eastern half, black sagebrush with mountain big sagebrush 
areas in swales and on north and east facing slopes. The 
climate is characterized by cold, snowy winters and warm 
summers. The predominant use of Parker Mountain is cattle 
and sheep summer grazing. Numerous antelope, sage grouse, 
Utah prairie dogs, and other wildlife are found in the area. 

METHODS 

Aerial Photo Interpretation 


The primary medium for preparing a 1:24,000 scale 
topographic map overlay of vegetation was high altitude CIR 
photography flown on Juiy 1-2, 1975. Film positive trans- 
parencies at 1:31,680 nominal scale were utilized. Mapping 
units were identified by examining the following: the 

color, texture, and patterns on the photographs; hydrologic 
features; topography (from topographic maps and stereoscopic 
viewing as needed); and ecological context. Scale adjust- 
ments and photographic displacement corrections were 
accomplished primarily by reference to U.S.G.S. orthophoto 
quadrangles, in addition to the use of a K&E Kargl carto- 
graphic projector. Interpretations were augmented by the 
use of a map from Landsat digital data (described below) . 



Three short trips to the study area in the summer and fall 
of 1982 were considered adequate for calibrating photo and 
Landsat interpretations with ground characteristics. 

Landsat Digital Data Analysis 

The Landsat multispectral scanner ("MSS") records light 
reflectance values for four spectral bands: green, red, and 

two bands of near Infrared light. MSS data represents light 
reflecting characteristics for the combined land cover ana 
terrain features within each picture element or "pixel," 
which covers approximately 1.1 acres of ground area. See 
U.S.6.S. (1979) for additional information regarding MSS data. 
Landsat data used in this study were recorded July 28, 1979. 

Landsat MSS data are analyzed statistically to detect 
light reflectance patterns which are sufficiently unique to 
make different ground cover types of interest consistently 
distinguishable (Hutchinson 1982) . The analytical approach 
used in this study, often referred to as an unsupervised 
classification method, began by examining the recorded MSS 
reflectance values for each pixel in the entire study area. 

From this search of individual pixels, statistics were 
generated which characterize pixel groups with similar spectral 
features. Next, a maximum likelihood classification routine 
was used to associate each pixel in the study area with one 
of the 46 spectral groups generated. 

The analysis next focused on detecting similarities and 
differences between spectral groups. A simple means of 
evaluating spectral characteristics is to plot each spectral 
group's mean reflectance value for the four MSS bands to form 
a diagnostic curve or "spectral signature." Since evaluating 
spectral signatures is often quite subjective, a more objec- 
tive technique was also applied. First, a principal compon- 
ents analysis of the mean values for each signature's four 
MSS bands reduced such data to factor scores for two compon- 
ents. Next, the factor scores were used in a cluster analy- 
sis which grouped spectral signatures according to a 
similarity index. Finally, the factor scores and group 
clusters were used in a discriminant analysis of the 
signatures. The two-dimensional scatter plot produced in 
the discriminant analysis allows one to receive a graphical 
view of signature relationships. (See Merola, et al. 1983 
for an example.) The use of discriminant analysis, based on 
MSS principal components and cluster analyses , in combination 
with examination of spectral signature plots and field experi- 
ence has been a key element in achieving good results from the 
unsupervised classification approach to Landsat data analysis. 

An additional and most vital dimension to the process 
of digital data analysis is calibrating spectral signatures 
with "ground truth." This is accomplished by assigning 
print symbols to each signature or signature group and print- 
ing maps which may then be registered to standard base maps 
or referenced to photographs and field study sites. In this 
study, a digital print map overlay was prepared to match the 
U.S.G.S. 7%-minute quadrangles (scale 1:24,000) mosaic of 
the study area. Calibration of spectral signatures with 
actual land cover types was accomplished primarily by use 




of th« v*g*t*eion map prepared from photo Interpretation, 
high altitude CIR photography, and field observations. The 
above-described process of Interpreting and combining 

S ectral signatures based upon signature curve similarity, 
acrlalnant analysis of the signatures and calibration of 
signature print symbols with photograph and ground observa- 
tions Is outlined In Rldd, et al.(1983J. 

The correspondence between Landsat spectral signatures 
and unique ground cover characteristics may be weak in some 
instances (Todd, et al. 1980). Landsat and pixel-analyzing 
computer algorithms perform robot- Ilk* functions and It Is 
often necessary to Introduce ancillary Information to 
Improve ground cover maps which are based solely on MSS data 
(Tom and Miller 1980). In this study, for example, the 
digital map calibration procass Indicated substantial spec- 
tral similarity, and therefore confusion, between sites 
dominated by a relatively lew growth form of big sagebrush 
and black sagebrush. This spectral similarity Is not 
surprising considering that both shrubs occupy similar 
ecological sites: generally on south and west facing slopes 

which are rocky and relatively drv. Such differences are 
also not evident on CIR photography, but must be ascertained 
In the field. However, black sagebrush appears to occupy 
this ecological site only on the western rim of the plateau 
and In areas to the east of a generalized 8,250 feet eleva- 
tion contour. To Improve the digital classification, the 
zone occupied by the short growth form of big sagebrush was 
digitized and an algorithm constructed to allow the detec- 
tion of differences between sagebrush species. Basically, 
the algorithm assigned each pixel with spectral signatures 
common to both species to different classes depending upon 
the location of the pixel with respect to the digitized zone. 

Other areas of spectral similarity were also addressed 
by the introduction of ancillary data. The surface geology 
of eastern half of the study area is predominantly older 
volcanic material in the north, and very recent volcanic 
flows in the south. The topography in the north Is charact- 
erized by a series of smooth ridges running in a southeast 
direction, whereas the southern area has more of a plateau 
character with various exposures. The combination of surface 
geology and topography differences between the areas has 
resulted in the confusion of big sagebrush, which grows in 
swales and northeast exposures in the north, with the black 
sagebrush signature of the south; since the black sagebrush 
areas which occur on southwest slopes in the north are 
spectrally different from black sagebrush on the recent 
volcanic flows in the south, there is little confusion bet- 
ween this spectral class and big sagebrush classes . The 
recent flows were digitized as separate units within the 
study area and new Landsat spectral class numbers were 
assigned to the signatures causing the confusion. 

Spectral similarity was also encountered in areas which 
are primarily bottomland loamy soils with mountain silver 
sagebrush or wetland vegetation cover. The majority of 
bottomland soils were digitized from an available S.C.S. 
soils map and spectral signatures not normally associated 



with Mountain silver sagebrush and wetlands were reassigned 
class nuabers to avoid confusion. 

RESULTS AND DISCUSSION 
CIR Photo Interpretation 

High altitude CIR photography proved to be an ideal 
photographic medium for the task of flapping rangeland 
resources; it provides high resolution prints with sore 
information and less displacenent than low altitude photo- 
graph a, and is relatively unaffected by atmospheric base 
which significantly scatters blue light. In addition, con- 
trasts between different vegetation types such as aspen and 
sagebrush are extremely vivid, whereas black and white photo- 
graphy often obscures such boundaries. CIR photography also 
generally produces greater discrimination between vegetation 
types than natural color photography because infrared light 
reflectance is highly sensitive to plant leaf shape and cell 
differences, as well as plant vigor. 

Despite the advantages offered by the CIR photography 
in this study, field observations and ecological interpreta- 
tions were vital in completing the mapping process. Generally, 
different ground cover types were found to be associated with 
distinct patterns of color, tone, and texture on the CIR 
photographs. A few circumstances led to confusion in Inter- 
preting the photos. For example, reddish-brown rocks on 
recent basalt flows often form a dominant feature which tends 
to produce similar CIR photo color- tone patterns, making 
black sagebrush areas indistinguishable from areas dominated 
by big sagebrush. Field observations and the nap from Landsat 
digital data helped to sort out the confusion. 

In a number of Instances, the ability to reference the 
Landsat digital print map significantly aided in the task of 
mapping rangeland cover from CIR photographs; the digital 
print map often flagged areas which sight otherwise have 
gone unnoticed because of subtle visual differences. Avail- 
able geologic, soils, precipitation, and topographic maps 
were also quite helpful. The rangeland cover types deline- 
ated from these procedures follow: black sagebrush; moun- 

tain big sagebrush, short growth form (dry, rocky sites); 
mountain big sagebrush, tall growth form (mesic sites); 
mountain silver sagebrush; wetland (rush/ sedge) ; aspen; 

Douglas fir; pinyon- juniper . 

Landsat Digital Data Analysis 

The methods applied in the analysis of Landsat data 
Initially expanded the number of spectral signatures before 
reducing the number of classes mapped to 20. As noted above, 
a total of 46 signatures were developed from statistically 
searching the study area for representative signatures. 
Partitioning the study area based upon elevation, geology, 
and soils, as described previously, led to the creation of 
more than 10 additional classes. The final selection of 20 
classes of rangeland cover represents a compromise between 
the goals of map simplification and preservation of meaningful 
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(or potentially meaningful) detail. Further feedback from 
digital print nap user* will determine whether the masher 
of cleaaee mapped ahould be expanded or reduced. It ie 
estimated that the Landsat map accuracy was increased by 
25X through the efforts to partition the study area based 
upon ancillary information. 

Landsat map accuracy wae assessed by randomly placing 
a grid, with vertices at ten-pixel row and column intervals, 
over the Landsat map and photo interpreting the nearest 
group of 4-10 pixels of a given class. A total of 830 pixel 
groups (average group sice was approximately 6) were 
examined, which represents a sample else of approximately 
12X of all pixels. Table 1 presents an error matrix for 
six levels of vegetation cover interpretation. Overall map 
accuracy is 89X, with the greatest amount of confusion 
associated with short growth form big sagebrush. Map 
verification by examining pixel groups probably produces a 
positive bias over verification of individual pixels because 
it leads to the checking of areas on photos which arc rela- 
tively homogeneous spectrally. However, simply looking at 
single pixels probably produces an opposite bias as a result 
of difficulty in achieving close registration between the 
Landsat up and photos. Sampling small groups of pixels is 
believed to be a good compromise, especially since pixels 
tend to occur as groups rather than as scattered individuals. 
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Table 2 illustrates that the Landsat map, based on 
available field information and photo interpretation, contains 
fair to good accuracy for most of the 20 spectral classes. 

The overall accuracy is 74X, with 21X average omission error 
and 23X average comission error. Table 2 also provides a 
brief description of the rangeland cover associated with the 
20 classes. 

Integration of Pholo Interpretation and Landsat Maps 

The ave ability of two different maps of rangeland 
cover could x '.d map users to feel obligated to select the 
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■or* accurate of ch* two approach**. However, ch*r* la no 
n**d to oak* auch a choice when the fundamental dlffarence* 
between both napping procedure* are examined. Both napping 
technique* have trade-off* in tern* of apatlal and inter- 
pretive accuracies which makes direct cooparlaon of naps 
difficult: relative map accuracies must be Judged bp refer- 

ence to available ground truth and in light of the particular 
spatial and Interpretive accuracy specification* of each 
mapping project. 
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Photo interpretation forces the sapper to generalise 
•pat tally to avoid creating sap polygons which are too 
numerous and/or too snail. Selecting ■<«<■»■ sapping unit 
slses and slapllfylng the sap legend are necessary In pre- 
paring a visually lntarpretable sap where patterns say be 
detected. Landsat sapping Includes nuserous cells, which 
are equivalent to. but generally such snaller than, the 
photo interpreter's line-drawn polygons, thus offering the 
potential for Increased spatial sapping detail. However, 
information obtained for each pixel by the Landsat scanner 
is already spatially generalised (ca. one acre pixel else), 
which offsets this advantage sosewhat. Of course, the dig- 
ital sap forsat offers the advantages of autosatsd area 
calculations, editing, updating, etc. 

Interpretive generalising and error occurs with photo 
interpretation since vegetation boundaries are not always 
distinct but, as a practical natter, lines nust be drawn to 
complete nap polygons. Generally speaking, the areal 
polygons sapped say be considered to have relatively homo- 
geneous land cover. However, subtle but significant 
vegetation nixing nay occur within sons sap polygons which 
is either undetected by the interpreter or overlooked to 
avoid additional lines and labeling of polygons which would 
clutter up the sap. Landsat data is quantitative, and 
analysis of spectral signatures is generally sore objective 
than photo Interpretation. The Landsat spectral bands 
cover a narrower spectral range than most photographic 
emulsions and permit the analysis of single or multiple 
bands: aspects which often serve to simplify the process 

of associating light reflactance with ground cover. Error 
and generalicing in interpretation of classes occurs when 
the combinations of physical factors, which determine multi- 
spectral reflectance for different land cover types, produce 
similar spectral responses. Unless such situations can be 
corrected by the use of ancillary data to digitlca boundaries 
which avoid the confusion, mlsclasslflcatlon of pixels will 
occur . 

Consequently, although Landsat analysis and CIR photo 
interpretations are both corns of remote sensing, comparing 
tha products produced by both methods is similar to Judging 
the difference between an apple and an orange. In addition, 
it can generally be assumed that ground truth is rarsly 
available in such abundance as to permit good comparisons 
of the mapping approaches; Landsat maps are typically eval- 
uated largely by reference to air photo interpretations, 
which leads to errors through misregistration or misinter- 
pretation and bias that nay inflate the apparent accuracy 
of maps from photo interpretation. 

The vegetation map produced in this study from photo 
interpretation offers the advantages of allowing a user to 
easily detect general patterns and locate most vegetation 
boundaries with a high degree of accuracy. The main weak- 
ness of the photo Interpreted map is that only six major 
vegetation classes vers able to be mapped, and that environ- 
mental variations within map polygons are often generalised. 
The Landsat map overlay in large part compensates for these 
weaknesses in the hand-drawn map with its 20 cover classes 
and nearly one acre polygon size. 
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Identlflcstlon of Ecolotlcel Reims Sites 

The delineation of Che vegetation cover types noted 
above provides an indication of present forage composition, 
but, more importantly, it provides a primary means of 
identifying range sites. The plant cownlty occupying a 
given site is a "synthomater” of the total environment of 
that site. The biotic and abiotic components have, over 
time, led to the dominance of the existing vegetation. 

Since Parker Mountain has not recently had any widespread 
raj or disturbances, it may be assumed that the vegetation 
types mapped from CIS. photography and Landsat data are the 
best single indicator of distinct sites for which manage- 
ment prescriptions mar be developed. This is especially 
the case where available information regarding spatial 
environmental variables such as soils, geologv, precipita- 
tion, and topography have been Integrated in botli land 
cover mapping approaches. 

Kangs sltea are mapped from an analysis of physical 
indicators which ar* associated with the range site. In 
this study, ecological conditions encountered on Parker 
Mountain suggest that twelve major range sites are present 
in the study area. These ecologlcsl sites have not been 
supped per it, but may be identified for any portion of the 
study area by reference to Table 3, and available maps of 
vegetation, soils, and geology. 
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CONCLUSION 

This project has permitted a cloaa evaluation of the 
relative Merits of mapping rangeland resources fro* CUL photo 
Interpretation and fro* MSS digital data. Beat results are 
obtained when both approaches are used in tandea; each ap- 

G oach has certain inherent disadvantages which are to a 
rge extent corrected by utilising the other approach. 

Combined with soils and geology naps , the vegetation 
aaps permit an accurate means of Identifying ecological 
range sites. At a cost of approximately 90.15/acre, this 
approach produces significant improvements in accuracy and 
efficiency over labor 'intensive alternatives. The maps and 
other information genarated as part of this study are pres* 
entiy being used by Parker Mountain range ranagers to select 
range condition/ trend monitoring sites and plan range 
improvements. 
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ABSTRACT 

Contextual Analysis of tandsat Spectral Signatures 

By 

John A. Mercia 

Center for Remote Sensing and Cartography 
University of Utah Research Institute 
420 Chlpeta Way, Suite 190 
Salt Lake City, Utah 84108 


A sequence of three multivariate statistical routines are presented 
to aid In detecting similarities and differences in spectral signatures. 
These signatures are derived from the picture elements (pixels) of 
Landsst MSS digital data. Several applications of this technique will be 
presented. 

The technique presented In this paper assists In producing, from MSS 
data, a more precise computer classified land cover map. Using and under- 
standing added spectral classes leads to Improved precision of the 
classified land cover map. The sequence of statistical routines used are: 
principle components analysis, cluster analysis, and discriminant analysis. 
The signature, or mean light reflectance values for each spectral class. 

Is used as the Initial Input data. 

Principle components analysis Is used to transform the data set 
(signatures) Into Independent variables to remove collinearlty. The 
variables now meet the assumptions of the statistical routines that follow. 
The output data are factor scores for two components representing the 
original data set. The factor scores are then used In a cluster analysis 
to group spectral’ signatures according to a similarity Index. Finally, 
the factor scores and group clusters are used In a discriminant analysis. 

A two-dimensional scatter plot produced In the discriminant analysis allows 
one to obtain a graphical view of the roectral context In which a particular 
signature Is found. 

The use of discriminant analysis, In combination with examination of 
spectral signature plots and ground Information, has been a key element In 
achieving good results from the unsupervised approach to Landsat data 
analysis. 
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DETECTING AGRICULTURAL TO URBAN LAND USE CHANGE 
FROM MULTI -TEMPORAL MSS DIGITAL DATA* 

Merrill K. Ridd, John A. Merole, end Richard A. Jaynes 
Center for Remote Sensing and Cartography 
University of Utah Research Institute 
420 Chipeta Way, Suite 190 
Salt Lake City, Utah 84108 

ABSTRACT 

Conversion of agricultural land to a. variety of urban 
uses is a major problem along the Wasatch Front, Utah. Loss 
of Agricultural production is one problem; another is the 
changing water consumption patterns that result from the 
land use conversion. Although Landsat MSS data is a rela- 
tively coarse tool for discriminating categoriea of change 
in urban-size plots, its availability sometimes to the 
exclusion of other diagnostic data prompts a thorough test 
of its power to detect change. 

This paper presents the procedures being applied to a 
test area in Salt Lake County, Utah, where the land conver- 
sion problem is acute. Results will be presented during 
the conference session. The objectives are to determine 
not only change but the identity of uses before and after 
conversion ana to compare digital procedures for doing so. 
Class selections are predicated on water consumption 
categories. The question is, how well can MSS data be made 
to identify these class types, and at what minimum plot 
size? Several algorithms are being compared, utilizing 
both raw data and preproeessad data. Verification of 
results Involves high quality color Infrared (CIR) photo* 
graphy and field observation. Selection of Landsat dates 
was determined by available CIR data, in this case involv- 
ing a two-year time span, from 1979 to 1981. The two data 
sets have been digitally registered, specific change 
categories are being internally identified in the software, 
results tabulated by computer, and "change maps" are to be 
printed at 1:24,000 scale. This digital analysla should 
help develop an automated procedure for wider area applica- 
tion and for subsequent updating without the necessity of 
repeated CIR photography. 

INTRODUCTION 

Salt Lake County, which represents Utah's major popula- 
tion center, continues tc experience rapid urban growth. 

The impacts of urbanization on land use patterns and natural 
resources in the county are of particular interest to both 
state and local policy makers. The effects of urban develop- 
ment on a dwindling agricultural land base and water 


♦This project is being supported by the National Aeronautics 
and "pace Administration (NASA Grant NAGW-95), in cooperation 
with the Utah Division of Water Resources. 


resources oust be assessed Co allow a rational basis for 
policy formulation. 

There is a need to have the capability to Monitor land 
uae changes on an annual basis. Since it is impractical to 
expect that aerial photography may be obtained at such 
frequency, the logical choice is Landsat to provide repeated 
remote sensing data. Landsat data offers the additional 
advantage of being in digital for*, thus obviating the need 
to perform the laborious process of repeated preparing and 
digitising of line maps. 

The obj active of this study nay be stated as follows: 
to utilise color infrared (CXR) photography and Maps to 
calibrate a digital land use map from Landsat data, for the 
purpose of establishing a data base which may be updated in 
future years without the necessity of having repeatad 
aerial photography coverage. Since this research project 
was underway at the time this paper was prapared, most of 
the results will be presented at the fall A.S.P. Conference. 
This paper focuses on the methods which have been and are 
being applied to accomplish the objectives noted above. 

STRATEGY 

The basic strategy is to establish 1979 as the base- 
line year, from which all subsequent change would be 
calibrated. For the baseline year, the most effective and 
accurate digital classification would be developed for the 
land cover types and sires Involved. The spatial /spectral 
data would remain on file for all future date comparisons. 
Updating might then be accomplished with simpler change 
detection algorithms. 

The procedure for this research is to (a) experiment 
with a variety of classification algorithms and select the 
combination yielding the best results, utlng this to map 
the 1979 baseline stage, and (b) to experiment with various 
change detection algorithms to determine the most effective 
update procedure. 


METHODS 

The first step was to identify a representative test 
area within the county. A test block approximately 15 ^>y 
15 miles was found to contain all the major urban and 
agricultural types of Interest. High quality CIR photo- 
graphy (1:30,000 scale), within one month of the Landsat 
digital data set, is being used to assist in identifying 
and calibrating cover features for that base year. 

The digital procedures described below are being used 
to detect changes occurring between the base year (1979) 
and 1981. CIR photography (1:60,000 scale), within one 
month of the Landsat data set, is being used in the 
analysis of the 1981 Landsat data. 

The analysis and mapping performed in this study are 
being carried out at the facilities of the Center for 
Remote Sensing and Cartography, using digital data 


obtained through NASA'* Landsat aatalllta system Croat tha 
EROS Data Cantor. For data processing, tha Cantor is 
utilising tha "ELAS" package of cootputar software routines, 
developed by NASA's Earth Resources Laboratory, which is 
operational on the University of Utah Research Institute's 
PRIME computer . 

Tha pr inary rationale for performing digital process- 
ing of MSS data begins with tha assumption that different 
types of ground cover have different patterns of reflection. 
It is also assunsd that these spectral patterns are 
sufficiently unique to make different ground cover types 
consistently distinguishable from one another using 
statistical classification techniques (Hutchinson 1982) . 

Standard Digital Data Analysis 

Ground cover characteristics may be analysed using 
digital image data from a variety of sources: Landsat 

multlspeetral scanner (MSS) ; Landsat thematic mapper (TM) ; 
airborne thematic mapper simulator (TMS) ; NOAA satellite 
sensor, advanced very high resolution radiometer (AVHRR); 
etc. The description of methods in this section applies to 
the analysis of digital image data from any such source. 

Any digital image data must first be reformatted to 
make them compatible with processing hardware. Next, a 
grey level map Is produced which encompasses the study area 
and soma of the surrounding area. The size of the grey 
level map is determined by the number of features present 
to be used as ground control. As water bodies typically 
provide the sharpest and most spectrally distinctive 
features for control, a near Infrared band is generally 
used to produce the map. (See Stage 1 of Figure 1.) These 
ground control points are used later to geographically 
reference the digital data to remove the effects of earth 
curvature, spin, etc. (See Stage 3 of Figure 1.) At this 
point of data processing, the raw data may be modified by 
other analytical steps such as filtering, principal compon- 
ents analysis, etc. Thereafter, a program called "SEARCH" 
is utilized to generate statistics which characterize 
pixel groups having similar spectral features across the 
bands.; - (See Stage 2 of Figure 1.) SEARCH ls'a routine 
which is used to provide training statistics for a program 
called "MAXL," which classifies individual pixels into a 
class based upon each pixel’s highest statistical probabil- 
ity of belonging to a given decs. The training statistics 
are derived from blocks or windows of data which correspond 
to areas on the ground containing the cover types of 
Interest. 

Once the study windows are selected, the SEARCH program 
examines each contlnguous six scan line (Landsat pixel 
matrix "row") by six element (pixel matrix "column") block; 
if the spectral data within the six bv six block are too 
heterogeneous , the program will switch to the use of a three 
by three block of pixels. The statistics generated by 
SEARCH Include mean pixel light reflectance values for each 
of the four bands, a covariance matrix, and a priori values. 








A ate of statistics Is gsnsrttsd by SEARCH representing 
various dasaas of light roflactanca pattarns found In eha 
study araa "searched. Tha mean light raflsctanca valuas 
for each spectral band are plotted to form a signature 
which characterizes each class. SEA&CH thus "trains" MAXL 
to recognize different ground cover patterns as It places 
individual pixels Into classes. A knowledge of the manner 
In which different land cover features create spectral 
signatures, combined with tha analysis of aerial photography 
ana field checking of digital classifications, allow the 
researcher to provide an Interpretation of Landsat-derlved 
spectral classes,. Characterizing the various signatures 
according to ground information. 


After light signatures are produced, further efforts 
are directed toward finding those signatures which would 
most likely reflect the major land cover types of Interest. 
Stage 2 of Figure 1 illustrates several of tha steps 
utilized In making detailed studies of signatures. A signa- 
ture plot permits a substantial amount of Interpretation; 
spectral signature shape and magnitude of reflectance are 
diagnostic of land cover types. Generally, similarly shaped 
signature curves Indicate similar cover types while upward 
or dmmward shifts of similar curves Indicate differences 
in topography, amount of ground cover, or the amount of 
mixing with other cover types, changing the overall 
"brightness." 


At CRSC, there has evolved a somewhat "standard" set 
of statistical routines to deal with an often unmanageable 
number of spectral classes. The sequence Involves principal 
components_anAlysl 9 , cluster analysis, and discriminant 
analysis, leading to a final classification. This sequence 
of routines may be applied to raw data or modified raw data 
(e.g., filtered, enhanced, etc.). The remainder of this 
section describes the application of these routines to 
Landsat MSS data, as performed for the Salt Lake County 
study area data. The analysis of signatures from other 
digital image data would be similar. 


Spectral signatures are then studied statistically to 
detect similarities and differences, which are s.ot always 
distinguishable from the signature plot. First, a principal 
components :.analysis .is of the mean values for each signature's 
four MSS bands' reduces such data to factor scores for two 
components; typically bands 4 and 5 are combined Into one 
component (visible light) , and bands 6 and 7 combine to form 
the second (Infrared light). Hext, the factor scores are 
used in a cluster analysis which groups spectral signatures 
according to a similarity index. Finally, the factor scores 
and group clusters are used in a discriminant analysis of 
the signatures. A two-dimensional scatter plot produced in 
the discriminant analysis allows one to obtain a graphical 
view of the spectral context In which a particular signa- 
ture is found. (See Figure 2.) The discriminant analysis 
scatter plot, with the two axes representing visible and 
infrared light components, may be divided into regions or 
groups of signatures that correspond to similar ground cover 
types. This process is a vital link in allowing an often 
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Figure 3 


A theoretical model of agriculture- to-urban land 
use conversion. 




uaaia>|ubli number of signatures to bs combined into 
groups of similar signatures. This procedure allows the 
researcher a great deal of flexibility In performing Landsat 
digital analysis; a large amber of signatures are available 
and one may concentrate on the signatures of particular 
Interest, while signatures of lesser Interest nay be grouped 
together or omitted. 

A vital dimension to the process of digital data 
analysis Is calibrating spectral signatures with ground 
Information. By assigning a print symbol to each signature 
or signature group, a print map siy os prepared and 
registered to standard base maps (e.g., U.S.6.S. 7fc minute 
quadrangles) and referrenced to photographs and field study 
sites. In this way, signatures are calibrated with actual 
land cover types. The use of discriminant analysis, based 
on MSS principal components and duster enalyses, In 
combination with examination of spectral signature plots 
and field experience (as outlined In Stages 2, 3, and 4 of 
Figure 1) , has been a key element In achieving good results 
from the unsupervised approach to Landsat data analysis 
(e.g., Merola, et al. 1983; Ridd, et al. 1983). 

Modified Raw MSS Data Analyses 

The procedures described above have been applied to 
both raw and modified Landsat data. 


Filtered Data . Reformatted Landsat data were spatially 
"filtered to increase homogeneity of pixel spectral values 
within a given cover type. The filter routine uses an 
unweighted moving average, or low pass filter. (See Jensen 
and Toll 1982.) In general, data filtering Involves adjust- 
ing raw spectral values for a given pixel to reflect the 
average spectral reflectance of a three by three pixel 
matrix around that pixel. This process tends to smooth out 
subtle spectral differences within an area which often 
results In the enhancement of edges betwoen areas that have 
significant spectral differences (Harallck 1979) . 

Principal Components . As part of the ELAS program pack- 
age, there is a program which derives principal components 
from raw data. The program then outputs, on a pixel- by- 
pixel basis, the component score for the plxsl based on the 
principal component being used. This program does not 
rotate the components. If It did rotate the components, 
one would be able to derive two components in a single 
analysis, one for the two visible bands and one for the two 
infrared bands; the two components would then be orthogonal 
or uncorrelated to one another (Johnston 1978) . Since this 
is not currently possible, we derived our components in two 
separate analyses, one to derive the visible component and 
one to derive the Infrared component (Walker and Shu 1982) . 
Although these comronents are not orthogonal to one another, 
they are very close to being so, and should provide 
information about the use or this technique. 


Other differencing routine* ere being examined a* well. 

A band 7/band 5 ratio was performed to eooblne the attribute* 
of those two bands to detect change. Another technique 
being evaluated i* a ratio of principal components. The 
ratio of visible/ infrared was done on the two dates separate- 
ly, then the two dates were differenced to form a differenced 
image. Another differencing technique utilises a Multi- 
Temporal Vegetation Index (KTVI) , presented by Derenyi, at 
al. (1962). A Vegetation Index (VI) was computed for each 
date as follows: 


VI- 
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then, the VI for one date is differenced from the other: 
KTVI - VI(T2)-VI(T1) 

Each of these techniques has shown some promise in various 
parts of the study area. Some are more effective for 
conversion of dry farm, others for rangeland, and others 
for irrigated land. 

Pastclasslflcatlon Improvement 

Following spectral classification, selected andllery 
terrain features are digitized to stratify cover types 
otherwise confused by spectral analysis. In addition, 
certain masks are created for "fixed" cover types, such as 
parks, golf courses, and public spaces that are spectrally 
confused with desired classes. These stratified and 
masked data files remain in the operating system for use 
In subsequent date change detection. 

Verification 


The accuracies of the several methods being explored 
are assessed using 247 ground verification plots. These 
plots have been digitized and can be compared to any map 
produces. The comparison is performed by a program in ELAS 
(ACTB) Accuracy of Classification Table. The table which 
is produced is an error matrix that shows class frequencies, 
percentages, percent correct, omission errors, and comission 
errors as a result of the comparison between the verification 
data and the classified data. 

A Model of Land Use Conversion 


The scatter plot in Figure 2 shews the 44 signatures 
derived from filtered data in the Salt Lake County test 
area. Signature groupings have been determined by analysis 
of the scatter plot, coupled with a signature plot from 
SEARCH, a classified print map, and field and photo 
observation. Each croup is associated with a particular 
land cover type, relevant to the land conversion process. 

Figure 3 presents a graphic illustration of a theoret- 
ical model of agriculture-to-urban land conversion in Salt 
Lake County. Irrigated farm lands occupy the upper left 




Thematic HiWtt Simulator Pitt Analysis 


Tha thematic mapper simulator (IMS) 1* a multispectral 
so aimer flown from an aircraft to simulate tha TM data that 
will ba available from Lands at 4. Tha goal in this analysis 
la to compare MSS data from aarllar Landset satallltaa and 
IKS data to ldantlfy tha possible role of TM and its us a in 
tha future. 


It has bean stated by several authors that an improved 
spatial and spectral sensor, such as tha IMS sensor, doss 
not necessarily improve classification accuracies whan 

K ocessed tha same as MSS data. The problem, as stated by 
arton (1982), Is as follows: "To illustrate the problem, 

consider the task of distinguishing between commercial and 
residential land-use classes in [spatially] high-resolution 
remotely sensed data. A straightforward par-pixel spectral 
classification of tha data can identify only the spectrally 
dissimilar ground cover classes such as pavement, lawn, 
trass, roof. These components are comon to both the com- 
mercial and residential land-use classes, making it 
difficult to distinguish between them on a par-pixel basis. 

It has been suggested by several authors that the solution 
to this problem is the use of contextual information (l.e., 
the local frequency distribution of components surrounding 
each pixel." A contextual raclassifier is being developed, 
at this writing, to address this problem. Spatial and 
spectral contextual information are applied in the procedure. 
The spectral contextual information is supplied by the 
statistical analysis of the signatures, as dtscribed 
previously. The two-dimensional scatter plot produced in 
the discriminant analysis is based on discriminant scores 
produced in the analysis. These discriminant scores supply 
the spectral contextual information. The spatial informa- 
tion is supplied by a moving window, of specified dimensions, 
and a similarity matrix, derived from the discriminant 
scores, to evsluate the center pixel of the moving window. 
The center pixel is then reclassified if the specified 
parsueters are met. This procedure will be tested on MSS 
as well as TMS data. 


Change Detection Data Analysis: Image Differencing 

Image differencing Involves subtracting the digital 
values or one image date from another, on a plxel-by-pixel 
basis. The first step in image differencing is to register 
each image to UTM coordinates using ground control points. 

The results of the differencing are transformed into positive 
values by adding a constant, and a histogram is produced 
from the differenced image. As Jansen and Toll (1982) have 
stated: "A critical element of the image differencing method 

is deciding where to place the threshold boundaries between 
change and no change pixels displayed in the hlatogram." 
Jensen and Toll found that image differencing with band 5 
produced the best definition of agricultural -to-urban change 
for fall data in tha Denver area. The present study has, 
tentatively, shown greater promise using band 7 data from 
mid- summer. 



region of the scatter plot. Prom this point, e field 
cleared for urban uae dramatically ehifte to the upper right 
region of the plot. Incipient residential (Ri) initially 
responds as barren land. New residential areas (Rn) have 
had some homes and streets constructed, which reduces 
brightness but have little effect on greenness. Initial 
landscaping established residential (Re) areas, and 
eventual maturing of residential vegetation (Rn) draw 
spactral characteristics toward the agriculture resion of 
the plot. 
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